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Abstract. We review recent work on the phase structure of QCD at very high baryon 
density. We introduce the phenomenon of color superconductivity and discuss how the 
quark masses and chemical potentials determine the structure of the superfluid quark 
phase. We comment on the possibility of kaon condensation at very high baryon density 
and study the competition between superfluid, density wave, and chiral crystal phases 
at intermediate density. 

1 Color Superconductivity 

In the interior of compact stars matter is compressed to densities several times 
larger than the density of ordinary matter. Unlike the situation in relativistic 
heavy ion collisions, these conditions are maintained for essentially infinite pe- 
riods of time. Also, compared to QCD scales, matter inside a compact star is 
quite cold. At low density quarks are confined, chiral symmetry is broken, and 
baryonic matter is described in terms of neutrons and protons as well as their 
excitations. At very large density, on the other hand, we expect that baryonic 
matter is described more effectively in terms of quarks rather than hadrons. As 
we shall see, these quarks can form new condensates and the phase structure of 
dense quark matter is quite rich. 

At very high density the natural degrees of freedom are quark excitations and 
holes in the vicinity of the Fermi surface. Since the Fermi momentum is large, 
asymptotic freedom implies that the interaction between quasi-particles is weak. 
In QCD, because of the presence of unscreened long range gauge forces, this is 
not quite true. Nevertheless, we believe that this fact does not essentially modify 
the argument. We know from the theory of superconductivity the Fermi surface 
is unstable in the presence of even an arbitrarily weak attractive interaction. 
At very large density, the attraction is provided by one-gluon exchange between 
quarks in a color anti-symmetric 3 state. High density quark matter is therefore 
expected to be a color superconductor 

Color superconductivity is described by a pair condensate of the form 

c(> = (ip T cr D \ C TFTp). (1) 

Here, C is the charge conjugation matrix, and .Tp, \ci t f are Dirac, color, and 
flavor matrices. Except in the case of only two colors, the order parameter can- 
not be a color singlet. Color superconductivity is therefore characterized by the 
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breakdown of color gauge invariance. As usual, this statement has to be inter- 
preted with care because local gauge invariance cannot really be broken. Never- 
theless, we can study gauge invariant consequences of a quark pair condensate, 
in particular the formation of a gap in the excitation spectrum. 

In addition to that, color superconductivity can lead to the breakdown of 
global symmetries. We shall see that in some cases there is a gauge invariant order 
parameter for the U(l) of baryon number. This corresponds to true superfluidity 
and the appearance of a massless phonon. We shall also find that for Nf > 2 color 
superconductivity leads to chiral symmetry breaking and that quark matter may 
support a kaon condensate. Finally, as we move to stronger coupling we find that 
other forms of order can compete with color superconductivity and that quark 
matter may exist in the form of chiral density waves or chiral crystals. 



2 Phase Structure in Weak Coupling 
2.1 QCD with two flavors 

In this section we shall discuss how to use weak coupling methods in order 
to explore the phases of dense quark matter. We begin with what is usually 
considered to be the simplest case, quark matter with two degenerate flavors, 
up and down. Renormalization group arguments suggest anc ^ explicit 

calculations show [^]||, that whenever possible quark pairs condense in an s- 
wave. This means that the spin wave function of the pair is anti-symmetric. 
Since the color wave function is also anti-symmetric, the Pauli principle requires 
the flavor wave function to be anti-symmetric, too. This essentially determines 
the structure of the order parameter [ p^JT]| 

This order parameter breaks the color SU(3) — ► SU(2) and leads to a gap for up 
and down quarks with two out of the three colors. Chiral and isospin symmetry 
remain unbroken. 

We can calculate the magnitude of the gap and the condensation energy using 
weak coupling methods. In weak coupling the gap is determined by ladder dia- 
grams with the one gluon exchange interaction. These diagrams can be summed 
using the gap equation 

AM = 4^ Jd qo J dcos, ( -1^^ (3 ) 

i + icos 9 \ A(q ) 



1- cos + ^/(2/^)7 ^ q 2 + A(qo) 2 ■ 

Here, A(p ) is the frequency dependent gap, g is the QCD coupling constant and 
G and F are the self energies of magnetic and electric gluons. This gap equation 
is very similar to the BCS gap equations that describe nuclear superfluids. The 
main difference is that because the gluon is massless, the gap equation contains 
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a collinear divergence for cos 9 — ► 1 . In a dense medium the collinear diver- 
gence is regularized by the gluon self energy. For q — > and to leading order in 
perturbation theory we have 

with m 2 — N f g 2 n 2 / (4tt 2 ) . In the electric part, m 2 D = 2m 2 is the familiar Debye 
screening mass. In the magnetic part, there is no screening of static modes, but 
non-static modes are modes are dynamically screened due to Landau damping. 
We can now perform the angular integral and find 

At ^ 9 2 f , , ( _^_ \ ^(go) 

18tt 2 7 Vbo-?o|/ Vffo+^o? 

with & = 2567r 4 (2/iV/) 5 / 2 <? _5 . This result shows why it was important to keep the 
frequency dependence of the gap. Because the collinear divergence is regulated 
by dynamic screening, the gap equation depends on po even if the frequency is 
small. We can also see that the gap scales as exp(— c/g). The collinear divergence 
leads to a gap equation with a double-log behavior. Qualitatively 
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from which we conclude that A ~ exp(— c/g). The approximation (J6j) is not 
sufficiently accurate to determine the correct value of the constant c. A more 
detailed analysis shows that the gap on the Fermi surface is given by 

A ~ 512n\2/N f ) 5 / 2 ng- 5 exp ("^) ■ (7) 

We should emphasize that, strictly speaking, this result contains only an estimate 
of the pre-exponent. It was recently argued that wave function renormalization 
and quasi-particle damping give O(l) contributions to the pre-exponent which 
substantially reduce the gap Jl6[ . 

For chemical potentials /i < 1 GeV, the coupling constant is not small and 
the applicability of perturbation theory is in doubt. If we ignore this problem 
and extrapolate the perturbative calculation to densities p ~ 5po we find gaps 
A ~ 100 MeV. This result may indeed be more reliable than the calculation on 
which it is based. In particular, we note that similar results have been obtained 
using realistic interactions which reproduce the chiral condensate at zero baryon 
density 

We can also determine the condensation energy. In weak coupling the grand 
potential can be calculated from [EJ 

« = \ J {"tr [S(q)S(q)} + tr [So\q)S(q)] } , (8) 
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where S(q) and S(q) are the Nambu-Gorkov propagator and proper self energy. 
Using the propagator in the superconducting phase we find 

<») 

where the factor 4 = 2Nf comes from the number of condensed species. Using 
A Q ~ 100 MeV and fi ~ 500 MeV we find e ~ -50MeV/fm 3 . This shows that 
the condensation energy is only a small 0(A 2 //i 2 ) correction to the total energy 
density of the quark phase. We note that the result for the condensation energy 
agrees with BCS theory. The same is true for the critical temperature which is 
given by T c = 0.56zAo- 

2.2 QCD with three flavors: Color-Flavor-Locking 

If quark matter is formed at densities several times nuclear matter density we 
expect the quark chemical potential to be larger than the strange quark mass. 
We therefore have to determine the structure of the superfluid order parameter 
for three quark flavors. We begin with the idealized situation of three degener- 
ate flavors. From the arguments given in the last section we expect the order 
parameter to be color and flavor anti-symmetric matrix of the form 

< = Wc-fttf)- (io) 

In order to determine the precise structure of this matrix we have to extremize 
grand canonical potential. We find |l7|p^| 

A? 3 = AWtf - (11) 

which describes the color- flavor locked phase proposed in |L9|. Both color and 
flavor symmetry are completely broken. There are eight combinations of color 
and flavor symmetries that generate unbroken global symmetries. The symmetry 
breaking pattern is 

SU(3) L x SU(3) R x U(l) v -» 517(3) v. (12) 

This is exactly the same symmetry breaking that QCD exhibits at low density. 
The spectrum of excitations in the color-flavor-locked (CFL) phase also looks 
remarkably like the spectrum of QCD at low density pof . The excitations can 
be classified according to their quantum numbers under the unbroken SU(3), 
and by their electric charge. The modified charge operator that generates a true 
symmetry of the CFL phase is given by a linear combination of the original charge 
operator Q em and the color hypercharge operator Q = diag(— 2/3, — 2/3, 1/3). 
Also, baryon number is only broken modulo 2/3, which means that one can 
still distinguish baryons from mesons. We find that the CFL phase contains an 
octet of Goldstone bosons associated with chiral symmetry breaking, an octet of 
vector mesons, an octet and a singlet of baryons, and a singlet Goldstone boson 
related to superfluidity. All of these states have integer charges. 
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With the exception of the U(l) Goldstone boson, these states exactly match 
the quantum numbers of the lowest lying multiplets in QCD at low density. In 
addition to that, the presence of the U(l) Goldstone boson can also be under- 
stood. The 17(1) order parameter is ((uds)(uds)). This order parameter has the 
quantum numbers of a + AA pair condensate. In Nf = 3 QCD, this is the most 
symmetric two nucleon channel, and a very likely candidate for superfluidity in 
nuclear matter at low to moderate density. We conclude that in QCD with three 
degenerate light flavors, there is no fundamental difference between the high 
and low density phases. This implies that a low density hyper-nuclear phase 
and the high density quark phase might be continuously connected, without an 
intervening phase transition. 



2.3 QCD with one flavor: Color-Spin-Locking 



In this section we shall study superfluid quark matter with only one quark flavor. 
As we will discuss in more detail in the next section, our results are relevant to 
quark matter with two or three flavors. This is the case when the Fermi surfaces 
of the individual flavors are too far apart in order to allow for pairing between 
different species. 

If two quarks with the same flavor are in a color anti-symmetric wave function 
then there combined spin and spatial wave function cannot be antisymmetric. 
This means that pairing between identical flavors has to involve total angular 
momentum one or greater. The simplest order parameters arc of the form 

= (e 3ab ip a C~ r ifj b ), $ 2 = (e 3ab ip a Cqip b ). (13) 



The corresponding gaps can be determined using the methods introduced in 



section gjj. We find A($ 1>2 ) = exp(-3ci, 2 ) A) with a = -1.5 and c 2 = -2 |§|J. 
While the natural scale of the s-wave gap is A = 100 MeV, the p-wave gap is 
expected to be less than 1 MeV. 



The spin one order parameter (13) is a color-spin matrix. This opens the 



possibility that color and spin degrees become entangled, similar to the color- 
flavor locked phase ( JTl| ) or the B-phase of liquid 3 He. The corresponding order 
parameter is 

* csl = 6?(e abc ij b C (cos(m + sin(/3) 7i ) (14) 

where the angle (3 determines the mixing between the two types of condensates 
shown in (|l3|). In reference we showed that the color-spin locked phase ( pr) ) 



is favored over the "polar" phase (13) 



In the color-spin locked phase color and rotational invariance are broken, but 
a diagonal SO(3) survives. As a consequence, the gap is isotropic. There are no 
gapless modes except if /3 takes on special values. The parameter /3 is sensitive 
to the quark mass and to higher order corrections. In the non-relativistic limit 
we find (3 = tt/4 whereas in the ultra-relativistic limit (3 = ir/2. 
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3 The Role of the Strange Quark Mass and the Electron 
Chemical Potential 



So far, we have only considered the case of two or three degenerate quark flavors. 
In the real world, the strange quark is significantly heavier than the up and 
down quarks. Also, in the interior of a neutron star, electrons are present and 
the chemical potentials for up and down quarks are not equal. 

The role of the strange quark mass in the high density phase was studied 
in |^l],^2). The main effect is a purely kinematic phenomenon that is easily 
explained. The Fermi surface for the strange quarks is shifted by 5pf = fJ> — 
((j 2 — m 2 ) 1 / 2 ~ to 2 (2/i) with respect to the Fermi surface of the light quarks. 
The condensate involves pairing between quarks of different flavors at opposite 
points on the Fermi surface. But if the Fermi surfaces are shifted, then the pairs 
do not have total momentum zero, and they cannot mix with pairs at others 
points on the Fermi surface. If the system is superfluid then the Fermi surface is 
smeared out over a range A. This means that pairing between strange and light 
quarks can take place as long as the mismatch between the Fermi momenta is 
smaller than the gap, 

2 

vn 

A> t, < 16 > 

This conclusion is supported by a more detailed analysis pl[p2| . Since flavor 
symmetry is broken, we allow the (ud) and (us) = (ds) components of the CFL 
condensate to be different. The Nf = 2 phase corresponds to (us) = (ds) = 0. 
We find that there is a first order phase transition from the CFL to the Nf = 2 
phase, and that the critical strange quark mass is in rough agreement with the 
estimate ([TBI). 

For densities p ~ (5 — 10) po the critical strange quark mass is close to physical 
mass of the strange quark. It is therefore hard to predict with certainty whether 
superconducting strange quark matter in the interior of a neutron star is in 
the color-flavor locked phase. Neutron star observations may help to answer the 
question. In the CFL phase all quarks have large gaps, whereas in the unlocked 
phase up and down quarks of the first two colors have large gaps, strange quarks 
have a small gap, and the up and down quarks of the remaining color have tiny 
gaps, or may not be gapped at all. 

The effects of a non-zero electron chemical potential was studied in |^3| . If the 
electron chemical potential exceeds the gap in the up-down sector then up and 
down quarks cannot pair. In this case, the up and down quarks pair separately, 
and with much smaller gaps, in the one flavor phase discussed in section |l4]. 
Between the phases with (ud) and (uu) , (dd) pairing there is a small window of 
electron chemical potentials where inhomogeneous superconductivity takes place 
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4 Kaon Condensation 



The low energy properties of dense quark matter are determined by collective 
modes. In the color-flavor locked phase these modes are the phonon and the 
pseudoscalar Goldstone bosons, the pions, kaons and etas. Some time ago, it 
was suggested that pions [ p5|j2^ , |27|] or kaons @,|9| might form Bose conden- 
sates in dense baryonic matter. The proposed critical densities were close to 
nuclear matter density in the case of pion condensation, and several times nu- 
clear matter density in the case of kaon condensation. Since mesonic modes 
persist in the color-flavor-locked phase of quark matter, we can now revisit the 
issue of Goldstone boson condensation |30|. In particular, we shall be able to 
use rigorous weak coupling methods in order to address the possibility of Bose 
condensation in dense matter. 

Our starting point is the effective lagrangian for the pseudoscalar Goldstone 
in dense matter |]3l| , |32[ 

C eff = ^Tr [d Ed E^ + v^Ed^] - c [det(M)Tr(M- 1 Z') + h.c] . (16) 

Here, E £ £{7(3) is the Goldstone boson field, is the velocity of the Gold- 
stone modes and M = diag(m„, nid, m s ) is the quark mass matrix. The effective 
description is valid for energies and momenta below the scale set by the gap, 
oj,q <C A. The low energy constants can be determined in weak coupling per- 
turbation theory. The result is vl = 1/3 and |§|3|,||,||j36) 



f2 21-81og(2)^ 



f = — — ov * (17) 
U 18 2tt 2 ' 1 ' 

A ~ (18) 



2^ 2 fl' 

We can now determine the masses of the Goldstone bosons 



m T ± = c(m u + m d )m s , m K± = cm d (m u + m s ). (19) 

This result shows that the kaon is lighter than the pion. This can be under- 
stood from the fact that, at high density, it is more appropriate to think of the 
interpolating field E as 

Eij ~ eikiejmn$L,k$L,li>R,m' t t>R,n (20) 

rather than the more familiar Eij ~ 4>L,iil>R,j fHf- Using (^0|) we observe that 
the negative pion field has the flavor structure dsus and therefore has mass 
proportional to (m u + rrid)m s ]32] |. Putting in numerical values we find that the 
kaon mass is very small, m K - ~ 5 MeV at fl = 500 MeV and m^- ~ 1 MeV at 
fi = 1000 MeV. 

There are two reasons why the pseudoscalar Goldstone bosons are anoma- 
lously light. First of all, the Goldstone boson masses in the color-flavor-locked 
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Fig. 1. Electron chemical potential (solid lines) and kaon mass (dashed lines) in the 
color-flavor-locked quark phase. The two curves for both quantities represent a simple 
estimate of the uncertainties due to the value of the strange quark mass and the scale 
setting procedure. 

phase are proportional to the quark masses squared rather than linear in the 
quark mass, as they are at zero density. This is due to an approximate Zi chiral 
symmetry in the color-flavor-locked phase pj[ . In addition to that, the Gold- 
stone boson masses are suppressed by a factor A/ p. This is a consequence of the 
fact that the Goldstone modes are collective excitations of particles and holes 
near the Fermi surface, whereas the quark mass term connects particles and 
anti-particles, far away from the Fermi surface. 

The fact that the meson spectrum is inverted, and that the kaon mass is 
exceptionally small opens the possibility that in dense quark matter electrons 
decay into kaons, and a kaon condensate is formed. Consider a kaon condensate 
(K~) = fxe" ¥i where p is the chemical potential for negative charge. The 
thermodynamic potential TL — pQ for this state is given by 



Here, p q = 3ps is the quark density, and x — p u /Pq and y — Pd/Pq are the up 
and down quark fractions. For simplicity, we have dropped higher order terms in 




(p 2 - m 2 K ) v 2 K + 0(v%) + pp q x 



1 



(21) 



12tt 
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the strange quark mass and neglected the electron mass. In order to determine 
the ground state we have to make J2l| ) stationary with respect to x,y,fi and 
vk- Minimization with respect to x and y enforces (3 equilibrium, while mini- 
mization with respect to /i ensures charge neutrality. Below the onset for kaon 
condensation we have vk = and there is no kaon contribution to the charge 
density. Neglecting m e and higher order corrections in m s we find 

In the absence of kaon condensation, the electron chemical potential will level off 
at the value of the electron mass for very high baryon density. The onset of kaon 
condensation is determined by the condition /i = mpc . At this point it becomes 
favorable to convert electrons into negatively charged kaons. 

Results for the electron chemical potential and the kaon mass as a function of 
the light quark Fermi momentum are shown in Fig. [l]. In order to assess some of 
the uncertainties we have varied the quark masses in the range m u — (3—5) MeV, 
md — (6 — 8) MeV, and m s = (120—150) MeV. We have used the one loop result 
for the running coupling constant at two different scales q = /i and q = fi/2. The 
scale parameter was set to Aqcd = 238 MeV, corresponding to a s (m T ) — 0.35. 
An important constraint is provided by the condition m s < \/2ppA discussed 
in the previous section. We have checked that this condition is always satisfied 
for pf > 500 MeV. Figure |l| shows that there is significant uncertainty in the 
relative magnitude of the chemical potential and the kaon mass. Nevertheless, the 
band of kaon mass predictions lies systematically below the predicted chemical 
potentials. We therefore conclude that kaon condensation appears likely even for 
moderate Fermi momenta pf — 500 MeV. For very large baryon density /i — > m e 
while m,K — > and kaon condensation seems inevitable. 

5 Chiral Waves and Chiral Crystals 

It is very important for the structure of compact stars to determine whether 
the quark matter core is in liquid or solid form. In the previous sections we 
discussed the case of weak coupling. In this case, particle-particle pairing is the 
only instability that needs to be considered. Nevertheless, in strong coupling, 
or if superconductivity is suppressed, other forms of pairing may take place. 
Obvious candidates are the formation of larger clusters or particle-hole pairing. 
Particle-hole pairing is characterized by an order parameter of the form 

(^(x)V(y)) = exp(zQ • (x + y))S(x - y), (23) 

where Q is an arbitrary vector. This state describes a chiral density wave. It was 
first suggested in |37j as the ground state of QCD at large chemical potential and 
large N c . This suggestion was based on the fact that particle-particle pairing, 
and superconductivity, is suppressed for large N c whereas particle-hole pairing is 
not. Particle-hole pairing, on the other hand, uses only a small part of the Fermi 



(22) 
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surface and does not take place in weak coupling. In the case of the onc-gluon 



exchange interaction these issues were studied in [38 . The main conclusion is 
that, in weak coupling, the chiral density wave instability requires very large 
iV c > 3. 

At moderate densities, and using realistic interactions, this is not necessarily 
the case. In particular, we know that at zero density the particle-anti-particle 
interaction is stronger, by a factor N c — 2, than the particle-particle interaction. 
In a Nambu-Jona-Lasinio type description this interaction exceeds the critical 
value required for chiral symmetry breaking to take place. For this reason we have 
recently studied the competition between the particle-particle and particle-hole 
instabilities using non-perturbative, instanton generated, forces [|39f . Our results 
are not only relevant to flavor symmetric quark matter at moderate densities, 
but also for the important case when there is a substantial difference between 
the chemical potentials for up and down quarks. As discussed in section || this 
disfavors ud-pairing, but it does not inhibit mm -1 and ddr 1 particle-hole pairing. 

Our results show that at low density the chiral density wave state is prac- 
tically degenerate with the BCS solution. Given the uncertainties that affect 
the calculation this implies that both states have to be considered as realistic 
possibilities for the behavior of quark matter near the phase transition. 

5.1 BCS Pairing 

In order to study competing instabilities we use the standard Nambu-Gorkov 
formalism, in which the propagator is written as a matrix in the space of all 
possible pair condensates. The BCS channel is described by the 2x2 matrix 

a _( (c*t4t> ( c ^ c -kl)\ = {G(k ,k,A) F(k ,k,A)\ 

CS "l(44t)(4 c -H)j -\F(ko,k,A)G(ka,-k,A)) ■ [M) 
The propagator has the form 

dsCS = - AA ( °-A at 1 ) ' (25) 

where 

Go = i l — ir> a ° = ITT \ -t ( 26 ) 

are the free particle propagator and its conjugate at finite chemical potential. 
Here, e k = d>k ~ /i q and S ek = sgn(efe)^ determines the pole position. From this 
equation we can read off the diagonal and off-diagonal components of the Gorkov 
propagator. The off-diagonal, anomalous, propagator is 

F ^ fe ' A ) = 71 I ■ r w,"l pn A 2 ■ ( 2? ) 

(ko - e k + lOe k ){ko + £fc + I0e h ) - A z 

The anomalous self energy A is determined by the gap equation 

-JL F( P0 ,p,A). (28) 
Here a pp is the effective coupling in the particle-particle channel. 
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5.2 Chiral Density Wave 



Using the same formalism we can also address pairing in the particle-hole chan- 
nel at finite total pair momentum Q. In the mean- field approximation, the full 
Greens function in the presence of a single density wave takes the form 



Govh — 




G(k ,k,Q,a) S(k Q ,k,Q,a) 
S(k , k, Q, cr) G(k , k + Q,Q, a) 



Go 1 



Again, we can read off the anomalous part of the Greens function 



S(k ,k,Q,cr) = 



—a 



(k - e k + iS ek )(k - e k+Q + iS ek+Q ) - a 



(29) 



(30) 



and the anomalous self energy a is determined by a self-consistency, or gap, 
equation 



a = (-i)a ph 



f d 4 p 
J (2tt) 4 



S(po,p,Q,°;tJ-q) ■ 



(31) 



Notice that the energy contour integration receives non-vanishing contributions 
only if 

e p e p+Q - a 2 < , (32) 

which ensures that the two poles in p arc in different (upper/lower) half-planes. 
This means that one particle (above the Fermi surface) and one hole (below the 
Fermi surface) participate in the interaction. 

The formation of a condensate carrying nonzero total momentum Q is as- 
sociated with nontrivial spatial structures. In the simplest case of particle-hole 
pairs with total momentum Q this is a density wave of wave length A = 2tt/Q. 
In three dimensions, however, we can have several density waves characterized 
by different momenta Q. In this case, the resulting spatial structure is a crystal. 
In general, the p-h pairing gap can be written as 



JnQj-r 



(33) 



where the Qj correspond to the (finite) number of fundamental waves, and the 
summation over |n| > 1 accounts for higher harmonics in the Fourier series. The 
matrix propagator formalism allows for the treatment of more than one density 
wave through a straightforward expansion of the basis states according to 



G = 



'k+Q x [l 
,t 



( c feT 



(cfe+Qxi 4t> ( C k+Q*l C l+Q x i) ( C fc+Qxi C k+Q y l) 

( C k+Q y l 4t> ( C k+Qyl 4+0*1 ) ( C k+Qyl C k+Q y l) 



(34) 
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Fig. 2. Dependence of the free energy (upper full line) and p-h pairing gap (dashed- 
dotted) on the number of waves ('patches') with fixed magnitude of the three- 
momentum \Qj\ — 0.8 GeV. The full line shows the value of the BCS ground state 
free energy. The results correspond to an instanton calculation with fi q = 0.4 GeV and 
N/V = 1 fm" 4 . 



The possibility of simultaneous BCS pairing can be incorporated by extending 
the Gorkov propagator to include both particle-hole and particle-particle com- 
ponents. In the following we will consider up to n w = 6 waves in three orthogonal 
directions with Q x = Q y = Q z and n = ±1, characterizing a cubic crystal. 

Note that in the propagators Go we do not include the contribution of anti- 
particles. This should be a reasonable approximation in the quark matter phase 
at sufficiently large /Lt g , when the standard particle- anti-particle chiral conden- 
sate has disappeared. At the same time, since our analysis is based on non- 
perturbative forces, the range of applicability is limited from above. Taken to- 
gether, we estimate the range of validity for our calculations to be roughly given 
by 0.4 GeV <; fi q <; 0.6 GeV. This coincides with the regime where, for the physi- 
cal current strange quark mass of m s ~ 0.14 GeV, the two-flavor superconductor 
might prevail over the color-flavor locked (CFL) state so that our restriction to 
Nf = 2 is supported. 

Solutions of the gap equations correspond to extrema (minima) in the energy 
density with respect to the gap a. However, solutions may exist for several values 
of the wave vector Q. To determine the minimum in this quantity, one has to 
take into account the explicit form of the free energy density. In the mean-field 
approximation, 

V 3 n(v q , Q,a)= f d 3 x ( ^ + (Q i (ia • V - 2a(x) q))) , (35) 




Fig. 3. Left panel: wave-vector dependence of the density wave free energy for one 
standing wave (full line: /2°" h , short-dashed line: O®^) in comparison to the BCS 
solution (long-dashed line: fluS s , dotted line: 0^ s ) at \i q — 0.4GeV. Right panel: 
wave-vector dependence of the density wave pairing gap (full line) compared to the 
BCS gap (long-dashed line). 



where V3 is the 3-volume. The first contribution removes the double counting 
from the fermionic contribution in the mean- field treatment. 

We have studied the coupled gap equations numerically. We do not find any 
solutions with simultaneous particle-particle and particle-hole condensates. This 
reduces the problem to the question whether the BCS or the density wave state 
is thermodynamically favored. The BCS solution A = 0.225 GeV is unique and 
has free energy of n B cs(n q = 0.4 GeV) = 2.3 • 1CT 3 GeV 4 . Here, we have 
neglected an irrelevant overall constant that does not affect the comparison with 
the density wave state. 

The situation is more complicated in the case of particle-hole pairing. Let 
us start with the 'canonical' case where the momentum of the chiral density 
wave is fixed at twice the Fermi momentum, Q = 2pp. In fig. || the resulting 
minimized free energy is displayed as a function of the number of included waves. 
The density wave solutions are not far above the BCS groundstate, with a slight 
energy gain for an increased number of waves. 

However, one can further economize the energy of the chiral density wave 
state by exploiting the freedom associated with the wave vector Q (or, equiv- 
alently, the periodicity of the lattice). For Q > 2pp the free energy rapidly 
increases. On the other hand, for Q < 2pp more favorable configurations are 
found. To correctly assess them one has to include the waves in pairs \k±Qj\ of 
standing waves (n w = 2, 4, 6, . . .) to ensure that the occupied states in the Fermi 
sea are saturated within the first Brillouin Zone. The lowest-lying state we could 
find at \i q — 0.4 GeV occurs for one standing wave with Q mm ~ 0.5 GeV and 
a ~ 0.21 GeV with a free energy fi ~ 2.3 • 10~ 3 GeV 4 , practically degenerate 
with the BCS solution. This density wave has a wavelength A ~ 2.5 fm. The 
minimum in the wave vector is in fact rather shallow, as seen from the explicit 
momentum dependence of the free energy displayed in fig. pi 
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